Quantum-control approach to realizing a Toffoli gate in circuit QED 



(N 

o 

a 

00 
(N 



43" 

-f— > ■ 



> 

od 

o 



X 



Vladimir M. Stojanovic, 1 A. Fedorov, 2 A. Wallraff, 2 and C. Bruder 1 

department of Physics, University of Basel, Klingelbergstrasse 82, CH-4-056 Basel, Switzerland 
2 Department of Physics, ETH Zurich, CH-8093 Zurich, Switzerland 
(Dated: January 31, 2012) 

We study the realization of a Toffoli gate with superconducting qubits in a circuit-QED setup 
using quantum-control methods. Starting with optimized piecewise-constant control fields acting 
on all qubits and typical strengths of Xy-type coupling between the qubits, we demonstrate that 
the optimal gate fidelities are affected only slightly by a "low-pass" filtering of these fields with the 
typical cutoff frequencies of microwave driving. Restricting ourselves to the range of control-field 
amplitudes for which the leakage to the non-computational states of a physical qubit is heavily 
suppressed, we theoretically predict that in the absence of decoherence and leakage, within 75 ns a 
Toffoli gate can be realized with intrinsic fidelities higher than 90%, while fidelities above 99% can 
be reached in about 140 ns. 

PACS numbers: 03.67.Lx, 03.67.Ac 



Superconducting (SC) qubits [l[ have come a long way 
since the realization that Josephson physics in SC circuits 
can be utilized to prepare well-defined few-level quantum 
systems [2|. Coupling such qubits is essential for quan- 
tum computation. The most successful approaches up 
to now rely on coupling all the qubits in an array to an 
"interaction bus," a central coupling element having the 
form of a transmission line with electromagnetic modes. 
In the circuit quantum electrodynamics (circuit QED) 
regime strong coupling between the qubits and the con- 
fined photons is realized ; the cavity photons induce 
a long-range coupling between the qubits. The combina- 
tion of transmon qubits @ and co-planar microwave cavi- 
ties represents the state-of-the-art of microwave quantum 
optics. 

The environmental degrees of freedom limit the time 
over which quantum coherence can be preserved. While 
they are similar to charge qubits, transmons have a much 
larger total capacitance such that the charging energy is 
significantly smaller than the Josephson energy. A small 
charge dispersion of the energy eigenstates leads to a re- 
duced sensitivity to charge noise and longer dephasing 
times (TV) tal. Recently, a significant progress has been 
achieved [7J,|£|, with T 2 times being increased by an order 
of magnitude from T 2 ~ 1 /is to T 2 ~ 20 (is. 

Given that two-qubit gates with SC qubits have been 
demonstrated with fidelities higher than 90% [§], a 
key challenge now is to realize three-qubit ones with 
shortest possible gate times. An example is the Tof- 
foli gate (controlled-controlled-NOT), which is relevant 
for quantum-error correction (lpj and has already been 
implemented with trapped ions [ll| and photonic sys- 
tems [IH with respective fidelities of 71% and 81%. Very 
recently, several groups realized a Toffoli gate using su- 
perconducting circuits 13-ll|. The Toffoli gate was im- 
plemented by a sequence of single- and two-qubit gates 
(direct approach). 

In this paper, we adopt an alternative approach and 
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FIG. 1: Lumped-element circuit diagram of three transmon 
qubits coupled to a superconducting transmission-line res- 
onator. The resonator serves as a coupling bus for the qubits 
and is also used for the readout of their states (red). Lo- 
cal flux lines (blue) allow for individual control of the qubit 
frequencies on the nanosecond time-scale. Microwave lines 
(green) are used to create control fields, each acting on its 
corresponding qubit. 



study the realization of a Toffoli gate with SC qubits in 
a circuit-QED setup (for an illustration, see Fig. [IJ by 
applying quantum-control methods [l6[ to the effective 
XY-type Hamiltonian of an interacting three-qubit ar- 
ray. We do so for realistic qubit-qubit coupling strengths 
and under typical experimental constraints on the qubit 
decoherence times. 

We first determine optimal piecewise-constant control 
fields and evaluate the resulting gate fidelities. Then we 
discuss how these fidelities are affected when the con- 
trol pulses are smoothened by eliminating their high- 
frequency Fourier components (spectral "low-pass" filter- 
ing). Our most important theoretical prediction is that 
within 75 ns a Toffoli gate can be realized with intrin- 
sic fidelities (in the absence of decoherence and leakage) 
higher than 90%, while fidelities higher than 99% can be 
obtained in approximately 140 ns. 

Methods of quantum control [l6| have been put for- 
ward in a number of theoretical proposals for realizing 
quantum logic gates with SC qubits 17 -Til, thus com- 
plementing studies that solely involve time-independent 
Hamiltonians [2^|. More fundamentally, recent studies 
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in operator (state independent) control 2l| have been 
focusing on interacting systems and employing the con- 
cept of local control [22], |23Jj . The essential idea is that 
systems such as coupled spin- 1/2 chains, models of inter- 
acting qubit arrays, can often be controlled by acting on 
a small subsystem. For instance, controlling only one end 
spin of an YYZ-Heisen berg chain ensures complete con- 
trollability of the chain 22, 23|. Yet, in view of the cur- 



aspects of the problem. Using the standard algorithm 



rently available few-qubit experimental setups [24( , more 
important than restricting control to only one qubit is to 
be able to carry out an optimal control-pulse sequence 
within times much shorter than T2. This is necessary in 
order to minimize the undesired decoherence effects 

Under the condition of resonant driving and assuming 
that the qubits are in resonance with one another, the 
effective (time-independent) YY-type (flip-flop) qubit- 
qubit interaction Hamiltonian is given by 



(see, for example, Ref. Il6l). it is straightforward to check 
that the dynamical Lie algebra of the system, generated 
by the skew-Hermitian operators —iHo, —i<Ji Xl and —ioiy 
(i = 1,2,3), has dimension 63 = d 2 — 1 (d = 8 is the di- 
mension of the Hilbert space of the system). This Lie 
algebra is isomorphic to su(d = 8) and the system is 
completely (operator) controllable. An arbitrary quan- 
tum gate can thus in principle be realized using properly 
designed control fields. 

Our goal is to find the time dependence of control fields 
(t) for realizing a quantum Toffoli gate. 
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We start our analysis with simple piecewise-constant con- 
trol fields [23| acting on all three qubits in alternation in 
the x- and y directions with control amplitudes f2 



x,n and 



n { y % {n=l,...,N t /2; 1 = 1,2,3). 

At t = control pulses are applied in the x direction 
to all three qubits with constant amplitudes il x \ during 
the time interval < t < T. The Hamiltonian of the 
system is then H x i = Hq + Y^h=i ^iViif Then y con- 
where a ix , <r iy , and a iz are the Pauli matrices. The sys- trol pulses with amplitudes are applied during the 
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i<3 



tern is acted upon by time-dependent Zeeman-like control 
fields described by the Hamiltonian 



H c (t) 
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[n^(t)a ix +^(t)a iy ] , (2) 



thus the system dynamics is governed by the total Hamil- 
tonian H(t) = H + H c (t). 

The control fields Cl x l \t) and Qy\t) can be imple- 
mented using arbitrary wave generators (recall Fig. [1}, 
which can produce an arbitrary signal with frequencies 
up to 500 MHz with minor distortions. In qubits based 
on weakly-anharmonic oscillators, leakage from the two- 
dimensional qubit Hilbert space (computational states) is 
the leading source of errors at short gate times [ill [2(| . 
This is especially pronounced if the control bandwidth 
is comparable to the anharmonicity 18]. We therefore 
impose an additional constraint on the control fields, viz. 
the condition that 



interval T < t < 2T, whereby the system dynamics is 
governed by H Vt \ = Hq + £^ =1 fl y l \(Jiy. This sequence 
of alternating x and y control pulses is repeated until N t 
pulses have been completed at the gate time t g = N t T. 
The time-evolution operator U(t = t g ) is then obtained 
as a product of the consecutive U x . n = cxp(—iH x ^ n T) 
and U Vt n = exp(—iHy tn T), where n = 1, . . . , N t /2. 

For varying choices of Nt and T, the 3N t control am- 
plitudes are determined so as to maximize the fidelity 
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tr[f/t(i ff )t/ TOFF ]| 



(4) 



<1„ 



t\fW^? + Wm 2 (3) 



is smaller than some threshold value for the transmon 
to be a well-defined two-level system. For typical an- 
harmonicities of transmon qubits (300 — 400 MHz) and 
values of f2 ma x of 100 — 130 MHz the error associated 
with this leakage should not exceed a few percent. While 
control schemes explicitly involving the higher levels of a 
physical qubit are in principle conceivable, in the present 
study we aim for simplicity and will in the following de- 
termine the control pulses for a system of qubits that are 
genuine two-level systems. 

Before evaluating numerically the optimal control 
pulses we would like to comment on the controllability 



where [/toff is the Toffoli gate. The numerical 
maximization over these amplitudes is carried out us- 
ing the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algo- 
rithm [53], a standard second-order quasi-Newton-type 
procedure. Based on an initial guess for the control 
amplitudes, the algorithm generates itcratively new se- 
quences of amplitudes such that in each iteration step the 
fidelity is increased, terminating when the desired accu- 
racy is reached. This procedure is repeated for multiple 
(~ 200) initial guesses to avoid getting trapped in local 
(instead of the global) maxima of F(t g ). The optimiza- 
tion is performed under the constraint f2 max < 130 MHz. 

While piecewise-constant control pulses are convenient 
as a starting point for a theoretical analysis, the actual 
pulse-shaping hardware cannot generate such fields with 
arbitrarily-high frequency components. Turning the time 
course of piecewise-constant control pulses into an opti- 
mized shape can be considered the central problem of 
numerical optimal control [28l |. 

We therefore perform spectral filtering of our opti- 
mal control fields. Quite generally, after acting with a 
frequency-filter function f(tu) on the Fourier transforms 
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T[ilf(t)] of the optimal fields nf\t) (j = x,y), one 
switches back to the time domain via inverse Fourier 
transformation to obtain the filtered fields flj(t): 

hf\t)=T- 1 [f(u J )Finf(t)}] (j = x,y). (5) 

In particular, we consider an ideal low-pass filter, which 
removes frequencies above the cut-off ujq and below — loq. 
In other words, /(w) = 9(lu + luq) —9{uj — ujq) , where 9{x) 
is the Heaviside function. When applied to our piecewise- 
constant control fields, the transformations in Eq. ([5]) can 
be carried out semi-analytically. They lead to 

, Nt/2 

^ " =1 (6) 

N t /2 K ' 

^(*) = - E tMt) - «2n-lW] , 

where a m (t) = Si[uio(mT — t)l (m G N) and Si(ar) = 
Jq (sint/t)dt stands for the sine integral. Based on 
Eq. ([5]) , we numerically determine the time-evolution op- 
erators corresponding to the filtered control fields using a 
product-formula approach (for details, see the Appendix 
in Ref. l23h . We then obtain the fidelities F(t g ) corre- 
sponding to the filtered fields from an analog of Eq. Q . 

Our numerical results are summarized in Table [I] It 
shows examples of calculated Toffoli-gate fidelities for 
optimal piecewise-constant control fields (F) and their 
low-pass filtered versions (F), for two different high- 
frequency cutoffs (u) = 500 MHz and u; = 450 MHz). 
The last column of the table shows the maximum fi max 

of ^j0\t)} 2 + [n { y\t)} 2 over all qubits and all times 
and obeys the constraint discussed after Eq. ©. The 
fidelities corresponding to the piecewise-constant fields 
are virtually unaffected by the filtering process, since the 
highest frequencies achievable by current pulse-shaping 
hardware are much larger than the qubit-qubit coupling 
strengths. 

To make contact with possible experiments, we now 
assume J = 30 MHz and J i2 = J23 = 6J13 = J [25j |. 
As can be inferred from the table, a Toffoli gate can be 
realized in 2.25J -1 = 75 ns with a fidelity higher than 
90%, while fidelities larger than 99% can be reached for 
gate times of around 4.18J -1 = 140 ns. Examples of 
optimal x and y control fields on all three qubits are 
shown in Fig. [3J 

For fixed total time t g = N t T higher fidelities are ob- 
tained for larger N t (i.e., smaller T), and the same is 
true of the robustness of these fidelities to random er- 
rors in the control-field amplitudes [23j. However, for 
T smaller than some (nonuniversal) threshold value, it 
becomes impossible to reach high fidelities (F > 90%) 
without violating the constraint f2 max < 130 MHz. 



TABLE I: Examples of calculated intrinsic Toffoli-gate fideli- 
ties for optimal piecewise-constant control fields (F) and their 
low-pass filtered versions (F), both corresponding to gate 
times t g given in units of J -1 where J = 30 MHz. The two dif- 
ferent values shown for F and r2 max correspond to respective 
high-frequency cutoffs of uio = 500 MHz and ujq — 450 MHz 
(in brackets). 















t g [ns] 


F[%] 


F[%] 


H max [MHz] 



14 
12 
10 

20 
18 
16 
18 
30 
28 
18 
20 
24 
22 
22 
24 
22 
22 



In reality, the fidelity loss resulting from decoherence is 
inextricably linked to the particular experimental setup 
and noise sources present in it. The errors due to de- 
coherence certainly depend sensitively on the total gate 
time, which is minimized in our approach by the interplay 
of always-on interactions between the qubits and time- 
dependent control pulses acting on all qubits. Quite gen- 
erally, the suppression of the gate fidelity due to decoher- 
ence is approximately given by the factor exp(— tg/T^), 
determined by the ratio of the gate time t g and the deco- 
herence time T2 0] . It is therefore quite encouraging that 
the required times we find for high-fidelity (F > 90%) 
realizations of the Toffoli gate (t g ~ 75 ns) represent a 
rather small fraction of the newly achieved decoherence 
times (T2 ^5 — 20 /j,s). Remarkably, this is better than 
achieved experimentally (with T2 ~ 1 fis) for two-qubit 
gates (t g ^ 30 — 60 ns). 

To summarize, employing methods of quantum opera- 
tor control we have investigated the feasibility of realizing 
a quantum Toffoli gate with superconducting qubits in a 
circuit QED setup. Our calculations indicate that within 
75 ns a Toffoli gate can be realized with intrinsic fideli- 
ties higher than 90%, while fidelities larger than 99% 
require gate times of about 140 ns. A particularly ap- 
pealing feature of our approach is that it does not make 
a principal difference between two- and three qubit gates, 
in contrast to the more conventional approaches in which 
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FIG. 2: Piecewise-constant and filtered (ojq = 500 MHz) con- 
trol fields acting on the three qubits for J = 30 MHz, i.e., a 
gate time t g of 4. 18 J -1 = 139 ns. 



a Toffoli gate is realized through several two-qubit and 
single-qubit gates. Our method can be used for realiz- 
ing other three-qubit gates (e.g., the Fredkin gate) and 
can also be straightforwardly generalized to an arbitrary 
number of qubits. 

Our study can be extended using more sophisticated 
iterative schemes for finding optimal control pulses, and 
could form the basis of an open-loop type approach tak- 
ing into account higher qubit levels and experimental un- 
certainties. The reduced gate times are likely to simplify 
the realization of three-qubit Toffoli gates and lead to a 

IMi- 



higher fidelity than the direct approach 
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